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^ . ABSTRACT 
\0 • 

I We investigate the linear polarization in the light of extrasolar planetary sys- 

^ I terns that may arise as a result of an occultation of the star by a transiting 

' planet. Such an occultation breaks any spherical symmetry over the projected 

^ ■ stellar disk and thus results in a non-vanishing linear polarization. This polar- 

00 ' ization will furthermore vary as the occultation progresses. We present both 

^ ' analytical and numerical results for the occultation of G-K-M-T dwarf stars by 

O ■ planets with sizes ranging from the one of Earth to two times the size of Jupiter. 

■ We find that the occultation polarization may result in an observable signal and 

I ' provide additional means to characterize various parameters of the system. A 

^ ■ particularly interesting result is that, for the later spectral types (i.e., smaller 

^ ■ stellar radii), this polarization signature may be observable even for Earth- like 

• planets. This suggests polarization as a possible tool to detect such planets. 

r~j ■ Departs from symmetry around mid-transit in the time dependence of the polar- 

j_| ■ ization signature may provide an estimate of the orbital eccentricity. 

Subject headings: polarization — occultations — Sun: general — planetary sys- 
tems 



1. Introduction 

In the wake of the first discovery of an extrasolar planet (ESP) orbiting a solar-type 
star by Mayor & Queloz (1995), over 150 ESPs have been found to date (see Ferraz-Mello 
et al. 2005, for a review). The vast majority of the discoveries have used the radial velocity 
method, in which high accuracy and precision Doppler measurements of the stellar spectrum 
yield the influence of the planet on the stellar motion. Among these ESPs, one (HD209458b) 
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has been observed to transit its star. As of May 2005 another 6 ESPs have been discovered 
photometrically through transits ^. Recently, Charbonneau et al. (2005) and Deming et al. 
(2005) reported the first direct measurement of the radiation of an ESP using the Spitzer 
Space Telescope. 

In addition to estimates of the ESP's period and, if the stellar mass is known, the semi- 
major axis of the planet orbit, the radial velocity curve provides a value for the minimum 
mass of the ESP, msini, and the eccentricity of the orbit. For a transiting planet, the 
inclination i can be estimated, and hence the ESP's mass, in addition to the planet's size 
and atmospheric contents (e.g., Alonso et al. 2004). 

This lopsided balance in favor of the radial velocity method as a means to detect ESPs 
may be turned around in the near future with the forthcoming satellite missions. As an 
example, the COROT mission is expected to detect hundreds of new planets using the pho- 
tometry method to measure occultations by transiting ESPs with micromagnitude accuracy. 
Such accuracy should also allow the detection of non-transiting ESPs through the added 
variation of their reflected hght along their orbit (Hough & Lucas 2003). 

The use of polarimctry as a tool to detect and study ESPs has been already expounded 
(Seager et al. 2000; Saar & Seager 2003; Hough & Lucas 2003; Stam et al. 2004). 
The main motivation in those works was precisely that the ESP will contribute to the hght 
beam with reflected starlight, at least in the optical and near infrared up to 3 /im (Stam 
et al. 2004). In the case of an unresolved star-ESP system, the observed polarization is 
expected to be rather small, ^ 10~^(= 10~^%) or less. This is because the unpolarized 
contribution from the star will be so much greater than the planet's reflected light (Seager 
et al. 2000; Saar & Seager 2003), the degree of polarization being essentially given by 
the ratio of reflected to total (= polarized + unpolarized) light. The direct detection of a 
resolved ESP via its reflected light has been detailed by Stam et al. (2004). Of course, for 
a planetary atmosphere the term scattering (by molecules and dust particles, for instance) 
of starlight might be in general more appropriate than reflection, at least for Jupiter-like or 
gaseous ESPs. 

Nevertheless, polarimetry is a relative measurement and high accuracy is possible. His- 
torically, the flrst mention of using highly accurate polarimetry for ESP detection seems to 
have been made by Kemp (oral communication) who built an equipment capable of ^ 10~^ 
accuracy, given enough photons (Kemp et al 1987). Hough & Lucas (2003) report on a 
polarimeter being built with a sensitivity better than 10~^ and aimed at detecting unresolved 
ESPs. At the same time, an equipment combining polarimetry, spectroscopy and adaptive 
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optics (CHEOPS, Feldt et al. 2003) has been proposed, which should achieve 10~^ fractional 
accuracy. All these are or will be ground-based equipment; a comparable photometric accu- 
racy can be typically obtained only from space. Light emitted by a spherical star will present 
a very small polarization. Kemp et al (1987) have measured a y-band linear polarization 
of less than 2 x 10^^ for the integrated solar disk, for example. However, when we resolve 
the solar disk, a polarization variation from center to limb will ensue (Leroy 2000). The 
center-to-limb variation in the polarization across any stellar disk originates from an exist- 
ing contribution of a scattering opacity to the total opacity in the atmosphere (Harrington 
1969). In hot stars, electron (Thompson) scattering is one such opacity. In the case of cool 
stars, Rayleigh scatter ing by atomic and molecular hydrogen and atomic helium provide 
the main source of scattering opacity. 

Barring any other source of asymmetry over the stellar disk (e.g., spots, flattening by 
rotation, etc.), the polarization pattern of a spherical star at a given limb angle will be 
tangential to the radial direction and the polarization value itself will typically increase from 
center to limb. This is because we tend to see more scattered, hence polarized, light towards 
the stellar limb. When a planet transits the star, it will then break this centro-symmetric 
pattern and cause a net linear polarization if we integrate the polarization across the stellar 
disk. 

In this paper we model and estimate the occultation polarization (henceforth OP) that 
results from the transit of an ESP. We consider stellar spectral types from G through T 
(brown dwarfs) and planet radii ranging from one Earth radius to two Jupiter radii. Wc first 
describe the problem and consider analytical approximations for the OP in Sect. 2. In Sect. 3 
we describe the simulation results for cool dwarfs and discuss the observational prospects as 
well as what can be learned from such observations. Finally, in Sect. 4 we summarize our 
results. 



2. Basic Equations 

Consider the geometry shown in Fig. 1-a, representing a planet orbiting a star with a 
circular orbit of radius a in the xy plane. At a given time, the position of the planet is given 
by its azimuth ip, measured from the x-axis. The direction to the observer is in the xz plane 
and has an inclination angle i with respect to the z-axis. The observer is at a distance d 
from the star (Fig. 1-b). 

Using this geometry we derive the expressions for the observed stellar flux and Stokes 
parameters Q and U. 
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2.1. Flux 

The stellar flux at O at a given wavelength is given by 

= I d(f) /o/(/i) cosasinada, (1) 
JO Jo 

where Iq is a constant, a* = arcsin [R^/d) is the angle at which a ray from O is tangent to 
the star, /x = cos 9 is the angle between the normal to the surface and the direction to the 
observer (Fig. 1-b) and /(/x) is the limb darkening law, deflned so that 

As will be evident later, a more convenient integration variable in Eq. (1) is p, the normalized 
impact parameter of a ray from O, given by 

p d sina 

i4 -ri* COS a 

In the limit a — > (or d >> R^,), p becomes 

d . 



P — 7^ sm a, 



R 



and equation (1) takes the form 



R' 



= d0 Pf{pi)dp. (2) 

From the relation p = \/l — jj?, also valid in the limit a — 0, we have = tt/q {Ri^/df , 
which is the familiar inverse square law for the stellar flux. 

Let us now calculate the flux absorbed by a planet with phase ip, observed from an 
inclination i. The projected coordinates of the planet center in the plane of the sky are, in 
the limit a ^ 0, 



, cos I ■ 

00 = arctan . (4) 



Pq = — y sin^ ip + cos^ ip cos^ i, (3) 
R* ^ 

fhn — arctan I 

tau'^ 

The flux absorbed by the planet is given by the integral of /o/(a*) over the circle projected 
by the planet on the stellar disk. From Eq. (2) we have 

?2 /•27r 



i^abs(Po,0o) = /o^y^ dcP'J^ P'f{p',cj>')dp', 
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where w = Rp/B^ and p' and (p' are polar coordinates of a system whose origin is (po, 0o)- 

A more useful quantity, which does not depend on d or Iq, is the observed normalized 
flux, defined as 

^* TT Jq Jq 



2.2. Polcirization 

We assume the limb polarization is a function only of P = P (/x), and the polarization 
direction is centrosymmetric with respect to the direction between the center of the star and 
the observer. Thus, for a given p and 0, the Stokes parameters are given by 

Q(p,0) = Jo/(/x)P(/x)cos(20), (6) 
C/(p,0) = -Jo/(/i)P(/i)sin(20). (7) 

It follows from Eq. (5), (6) and (7) that the observed normalized Stokes parameters, q — 
Q/Fahs and u = U/F^^s, are 

q{Po: 0o) = - / d0' / p'fip', <l>')P{p', 0') cos(20)dp', (8) 
71" Jq Jo 

1 /•27r pw 

u{Po: 0o) = — / / p7(p', <I>')P{P': 00 sin(20)dp', (9) 

Jo Jo 

Note that in the above expressions, the argument of the sine and cosine is 20, without the 
prime. Note, also, that the expressions above for the observed flux and Stokes parameters 
are not valid for the ingress and egress phases of the transit. 



2.3. Analytical approximations 

In this section we derive a simple analytical approximation for q and u, valid when 
w <^ 1. In this case, the argument of the integrals in Eqs. (8) and (9) is approximately 
constant 

/(p',0')^/(Po), 
P{p',<P')^P{po), 

and the observed normalized Stokes parameters take the following simple form 

QiPo, (Po) ~ w'^fipo)P{Po) cos(2(/)o), 
u{poAo) ~ -tyV(Po)-P(Po)sin(20o)- 
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The observed OP in this approximation is given by 

piPo) = [qiPof + uiPof] ^ w'fiPo)P{Po). (10) 

The polarization depends only on the impact parameter, because of the centro-symmetric 
pattern of the stellar limb polarization. 

Let us now obtain an expression for the maximum OP during the transit. Because the 
limb polarization usually grows towards the limb, we expect that the maximum polarization 
occurs near the interior contact of either the ingress or egress phases. The impact parameter 
of the planet center at this phase is po = 1 — w, hence the maximum polarization is 

Pms.^W^f{l-w)P{l-w)- (11) 

The above expression shows that the maximum polarization is controlled by two quan- 
tities with opposite effects. The first quantity is the ratio of the sohd angles subtended 

by the planet and the star (w^); the larger this ratio, the larger the flux absorbed by the 
planet and the larger the polarization. The second quantity is the limb polarization at the 
impact parameter of the planet center at the contact. Here, the larger the size of the planet, 
the lower the impact parameter and, thus, the lower the polarization. As we discuss in the 
following section, the existence of these two competing effects has important consequences 
on the maximum OP produced. 



3. Numerical Results 

3.1. Monte Carlo Code 

In the following sections, we present numerical polarization results for representative 
cases of G, K, M and T dwarfs. We used a Monte Carlo integrator for solving Eqs. (8) 
and (9). For a given configuration of the star-planet system (i.e., B^, Rp, a, and the 
procedure adopted is as follows. A photon packet (PP) is emitted in a random position on 
the stellar surface and in the direction towards the observer. The PPs are given the following 
weight 

e. = ^, (12) 

where is the stellar luminosity. The above expression means simply that the weight is 
proportional to the probability of the PP being emitted towards the observer. In addition 
to the weight, the PP is characterized by its linear polarization 

= QP(Ai)cos(20), (13) 
C/i = -QP(/x)sin(20). (14) 
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If the PP is emitted in the hemisphere opposite to the one facing the observer it is given 
weight zero, because it can not reach the observer. If, on the contrary, the PP is emitted in 
the hemisphere facing the observer, it is checked whether its trajectory crosses the planet. 
If it does, the PP will be absorbed by planet and, thus, is also given weight zero. 

This procedure is repeated times and, at the end of the simulation, the flux and 
Stokes parameters will be given by the following estimators 



F^ = J2w„ (15) 



N 

- — ^ — > (J-oj 

u, = (17) 

where the subscript ip indicates the dependence of the estimators with planet position. F, q 
and u are calculated for several planet positions in order to simulate the temporal progression 
of the transit. 

The above procedure is extremely simple to implement and it has the advantage of 
avoiding dealing with the complications that arise on the ingress and egress phases of the 
transit, for which only part of the planet disk blocks stellar light. In addition, other physical 
effects (such as stellar spots) can be easily implemented. 



3.2. Limb Darkening and Polarization 

For the hmb darkening law we use the formula proposed by Claret (2000) 

fifx) - 1 - ai(l - ii'/^) - a2(l - A^) - a3(l - /^'/') - ^4(1 - /x'), (18) 

where the are constants that depend on wavelength, spectral type, surface gravity and 
metaUicity. In Table 1 we hst the we used for each spectral class, as well as the radius 
and effective temperatures adopted. Solar metallicities are assumed. 

In contrast to the limb darkening, which is relatively well-known for normal stars, the 
limb polarization is well-known only for the Sun, for which it was extensively studied both 
observationally and theoretically (e.g., Fluri & Stenflo 1999; Leroy 2000; Faurobert et al. 
2001; Faurobert & Arnaud 2003). For other spectral types, however, to our knowledge only 
a handful of theoretical studies has been made (Harrington 1969, 1970; Magalhaes et al. 
1986; Collins 1988). 
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For the Sun, we adopt a limb polarization of the following form (Fluri & Stenflo 1999) 

where A; is a measure of the width of the polarization profile (i.e., the larger the k, the 
fastest the polarization drops from the value at the limb). Both Pi and k are a function 
of wavelength. Fluri & Stenflo (1999) studied the solar limb polarization between 4000 
and 8OOOA and found that the above expression matches reasonably well both models and 
observations. Values of Pi and k for three wavelengths are hsted in Table 2. 

For the later spectral types, whose limb polarization is not known, we use Eq. (19) and 
explore different values of Pi and k to study how the polarization signature varies with these 
parameters. 



3.3. G Dwarfs (the Sun) 

In the models shown in this and the following subsections we investigate planet radii 
ranging from 0.09 to 2 Jupiter radii, orbiting the star with a circular orbit of radius 0.04 AU, 
a typical distance for the so-called hot Jupiters (Silva & Cruz 2005). We point out that the 
only consequence on the following results, which arise from this choice of planet distance, is 
in the time-scale of the transit. 

The transit of ESPs with radii 1 it!jup and 2 it!jup in front of a Sun-hke star is depicted 
in Fig. 2 for three different inclination angles. The figure illustrates the size of the planet 
compared to the star as well as the path of the transit across the stellar disk for the adopted 
inclinations. 

The monochromatic normalized fiux and OP at 4600A as a function of time, resulting 
from the geometry of Fig. 2, are shown in Figs. 3 and 4, respectively. The dashed horizontal 
line in Fig. 4 indicates the best polarization sensitivity for current or planned broad-band 
polarimeters (fs 10~^, Hough & Lucas 2003). 

From Figs. 3 and 4 it is evident that the effects of inchnation angle are quite different 
for the fiux than for the polarization. For the fiux (Fig. 3), decreasing the inchnation angle, 
which shifts the transit path on the stellar disk to higher latitudes, makes the fiux curve 
shallower as a result of the limb darkening. Observationally, this makes photometric transit 
detections more difficult for smaller inclinations. 

In constrast. Fig. 4 shows that the shape of the OP changes dramatically with inchnation 
angle; it goes from a double-peaked profile for i — 90° and intermediate inclination angles 
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to a (nearly) single-peaked function for lower inclination angles. This happens because, as 
discussed in Sect. 2.3, the maximum OP occurs near the interior contacts of either the ingress 
and egress phases, for which the projected planet disk is close to the stellar limb. For the 
same reason, the maximum OP is not dependent on the inclination, provided the angle is 
not so small that only part of the planet disk transits in front of the star. Halfway in time 
accross a central {i = 90°) transit, the polarization is zero because the projected stellar disk 
becomes symmetrical, which does not happen for other inclination angles. 

The above has two important observational consequences. First, the fact the the maxi- 
mum OP is independent on the inclination means that low inclination angles will not hamper 
the OP detection, which may not be the case for photometric observations. Second, the 
strong dependence of the OP curve with inclination angle suggests polarization as a good 
indicator for this parameter. 

As far as polarization goes, the results of Fig. 4 represent the most unfavorable scenario 
we explore in this paper, because of the large size of the star compared to the size of the 
planet and the low values of the hmb polarization. Even so, the polarization levels in Fig. 4 
indicate that the OP ensuing from the transit of Jupiter-like or larger planets in front of G 
dwarfs should be measurable. Larger OPs for G Dwarfs are expected in the bluer regions 
of the spectrum, because the values of Pi increases fast with decreasing wavelength (Leroy 
2000). For instance, at A = 3000A Pi is about 10 times larger than at 4600A, and the OP 
in this spectral region would, thus, be proportionally larger. 

Large OPs are also expected at the wavelengths of the resonant lines of several elements, 
such as Ca and Sr. For instance, Bianda et al. (1999) measured the solar limb polarization for 
the resonant line Ca 1 4227A and found the limb polarization to be quite large (Pi = 16.5%). 
Using these data, we constructed the OP for this wavelength and show it in Fig. 5. It can 
be seen that it is indeed large and well above the detection threshold. 



3.4. K, M, and T Dwarfs 

In this subsection we investigate the OP for later spectral types, for which higher polar- 
ization levels are expected. This is due to both the larger ratio between the planet and stellar 
radii and the higher limb polarization levels expected for those spectral types (Harrington 
1969). Since we do not know the limb polarization in these cases, we define the quantity 




which can be regarded as a polarization efficiency of the star-planet system. Given the 
adopted detection threshold of 10~^ and a value of Pi, it follows that the OP should be 
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detectable if e > 10~^/Pi. Based on the observed solar values (Table 2), a typical minimum 
value for Pi is ~ 0.001, and e should then be typically greater than about 0.001 if an OP is 
to be measured. 

In Fig. 6 we first illustrate how the limb polarization curve of Eq. (19) varies with k. 
Large values of k, representative of the solar values (Table 2), result in a very sharp decrease 
of the limb polarization. For k = 50, for instance, the polarization drops to 10% of Pi about 
0.02 stellar radii away from the limb. For latc-type giants, the limb polarization curve is 
probably associated with low values of A; (A; a; 1, see caption to Fig. 6). 

In Fig. 7 we show the polarization efficiency e for the central transit of a 1 R^^p planet 
in front of a M2 dwarf, for which we adopted a radius of 0.5 Rq (Table 1). The figure shows 
that the parameter k has a dramatic effect on the polarization efficiency. For k between 
and 5 e is quite large, resulting in observable OP levels {eP\) even for low (a 1%) Pi. For 
A; = 50 the efficiency is much smaller, but still large enough to produce OPs greater than 
the detection threshold, for large planets. 

The explanation for the strong dependence of e with k is straightforward. As discussed 
in Sect. 2.3, the peak of the OP occurs near the inner contact of either the ingress or egress 
phases, and is controlled by the limb polarization at the impact parameter of the planet 
center at the contact. From Eqs. 11 and 19 we can write the maximum efficiency as 

emax ~ W'f{l - W) ^ l'"^^^ « W^f{l - w) ^ ~ (21) 

^ ^ + A;Vl-(l-w)2 ^l + kV2^ ^ ' 

It is evident from the above expression that the maximum efficiency drops with increasing 
k. 

In Fig. 8 we show how the maximum polarization efficiency varies with the planet 
radius. The figure illustrates well the two opposite effects that control the OP, as discussed 
in Sects. 2.3 and 3.3. Increasing the planet radius enlarges the flux absorbed by the planet 
and makes the polarization go up. This increase, however, is hampered by the fact that the 
limb polarization at the impact parameter of the planet center at the inner contact decreases 
with increasing planet radius. This is the reason why the curves in Fig. 8 level off somewhat 
as Pp increases. In any case, all the polarization values implied by the efficiencies shown in 
Fig. 8 should be detectable. 

The dotted lines in Fig. 8 represent the analytical approximation for the maximum 

OP, Eq. 11, with the appropriate values for k and the limb darkening coefficients. The 
differences between the analytical and numerical results are small, which indicates that Eq. 11 
is very useful for providing an estimate of the maximum OP for a given configuration. 

In Figs. 9 and 10 we show the results for a K5 dwarf, for which we adopted a stellar 
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radius of 0.72 Rq. The results for a T dwarf with radius 0.2 Rq are shown in Figs. 11 
and 12. We note that theoretical atmosphere models for brown dwarfs (e.g., Tsuji 2002) 
indicate the formation of dust grains in such cool atmospheres. Those grains will provide 
an additional scattering opacity to the Rayleigh opacity, suggesting large values for the limb 
polarization for those objects. Such limb polarization needs further detailed modeling. 

Our results clearly demonstrate that the OP is within the reach of current polarimeters. 
In the worst scenario, corresponding to the red lines of Figs. 7 to 12 {k = 50), the efficiency 
ranges between about 10^"'' (i?p = 0.09 -Rjup) to 0.003 (i?p = 2 -Rjup), which corresponds to 
polarization levels of 10^^ to 3 x 10^^ if the Pi = 1%, and ten times those values if Pi = 10%. 
The best scenarios correspond to /c = (black lines) . In these cases the efficiency is between 
3 X 10~^ and about 0.1. For a limb polarization of 10%, this would correspond to polarization 
levels of 3 X 10~^ up to 1%. The maximum OP values for the hmiting cases studied here are 
summarized in Table 3, where we adopted Pi = 10%. 

Table 3 demonstrates that high precision polarimetry can be a means not only to study 
planetary systems with large (Jupiter size) planets, but also to detect Nepture-like planets 
(Pp 0.4 Pjup), and possibly Earth-hke planets, depending on the shape and level of the 
limb polarization. 

We also point out that the OP curves, exemphfied by Figs. 3, 4, 7, 9 and 11, are 
symmetric with respect to time because of the circular orbit adopted. An elliptical orbit 
will result in an OP time dependence which will be non-symmetrical with respect to mid- 
transit, the details of which will depend on orbital parameters such as the orbital major-axis 
orientation and eccentricity. Thus, observations of the OP may potentially provide important 
additional orbital information which may be difficult or impossible to obtain from the light 
curve of the event alone. In particular, knowledge of the eccentricity values for a sample 
such as the one that might result from the COROT mission should provide further insight 
into the dynamical evolution of the cxoplanet systems thus observed. 

An important point to consider is the temporal resolution that may be attainable for high 
precision polarimetric measurements. For a given absolute error dp, one requires typically 
~ 2/(Tp photons. For a 9th magnitude star, for instance, about an hour of exposure time on 
a 4 meter telescope will be needed to reach (jp 10~^. As a result, the adopted limit of 10~^ 
can only be reached for long period transits (several hours to days) or for very bright stars. 
For short transits (r^ 1 hour) and/or fainter stars the detection threshold for measuring 
polarization variabihty across the event may be larger than 10~^. On the other hand, as we 
pointed out above, transits that involve cool stars and/or large planets may result in OP 
values well above 10~^. Also, experiments such as the COROT satellite will likely detect 
long- duration transits. 
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We want to briefly discuss another possible source of variable polarization that might be 
misinterpreted as a planetary transit. Stellar spots may mimic a planet transit by removing 
part of the expected signal from some regions of the stellar disk, which breaks the spherical 
symmetry over the stellar disk, producing a polarization signal that could be similar to 
that of a transiting planet. Despite these anticipated similarities, it may be possible to 
separate any long-lived-spot polarization signature from the OP with detailed modeling, 
which is outside the scope of this paper. For instance, the net polarization as the spot moves 
into the observed hemisphere will likely depend on the temperature contrast between the 
spot and the underlying photosphere, as well as the detailed limb polarization from each 
location. Additionally, the spot may be intrinsically polarized by magnetic mechanisms 
(Saar & Seager 2003). Finally, it might be possible to observationally distinguish the 
two phenomena, spots versus oculting planets, from detailed differences in their respective 
photometric (and possibly spectroscopic) variations. 

Finally, tidal distortion induced by the planet on the star will originate a polarization 
signal that will vary with a period identical to the orbital period with an expected level of 
?a 10~''-10~^ (Hough & Lucas 2003). This interesting effect is expected to be important 
only for massive planets and is typically lower than the polarization signal discussed in this 
paper. 

4. Conclusions 

In this paper we studied the occultation polarization (OP) that ensues from the transit 
of an extra-solar planet in front of cool dwarfs. Following the derivation of some simple but 
useful analytical approximations, we presented in Sect. 3 numerical results for the OP for 
various planet sizes and stellar spectral types ranging from G to T. 

For the spectral type G, we used limb darkening and polarization from the Sun to 
calculate the monochromatic OP at 4600A and for the resonant line Gal (4227A). We found 
that even for such a (relatively) large star and small values of the limb polarization (8.7 x 
10^'^), the OP is generally large enough (> 10^^ = 10~"^%) to be observed with modern high 
precision polarimeters being built (Hough & Lucas 2003). As discussed in the text, larger 
values are expected for the OP of G dwarfs both in the bluer parts of the spectrum and in 
the resonant lines of several elements, for which the limb polarization is larger than in the 
visual continuum. 

Because of the lack of knowledge of the hmb polarization for the later spectral types 
(K to T), we adopted the formula of Eq. 19 and explored several values of k, the width of 
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the limb polarization curve. We found that in the least favorable cases (i.e., those of large 
k) the detection of the OP of Jupiter-like planets, down to Neptune-like planets, should 
be possible. In the most favorable situations (small k), our results suggest the interesting 
prospect of being able to detect Earth-like planets orbiting late type dwarfs and brown dwarfs 
via ground-based polarimetry. 

We also studied the properties of the solution with respect to the inclination angle of 
the observer. We found that, while flux observations of the transit may suffer from weaker 
signatures at smaller inclination angles (i.e., transit path along high stellar latitutes), the 
same docs not happen for the polarization, because the maximum OP is independent of i. 
Also, we found that the shape of the OP is much more dependent on i than the shape of the 
flux curve, suggesting that polarimetry can help better constrain this parameter. 

This work has been supported by Sao Paulo State Funding Agency FAPESP (grants 
01/12589-1 and 04/07707-3). AMM acknowledges support by the Brazihan agency CNPq. 
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a) b) 




Fig. 1. — Transity geometry. The observer's direction is in the xz plane and has an angle i 
with respect to the z-a.xis. 




Fig. 2. — Geometry of the transit of a 1 Rj^p ESP (sohd circle) and a 2 i?jup ESP (dashed 
circle) in front a Sun-like star, for three different inclinations. For comparison, the geometry 
is equivalent to 0.5 it!jup and 1 it!jup ESPs transiting a M2 dwarf and to 0.2 and 0.4 it!jup 
ESPs transiting a brown dwarf. 



- 16 - 




0.94 t , , , I , , , I , , , I , , , ^ 

-4-2 2 4 

t (h) 

Fig. 3. — Normalized flux at 4600 A for the transit of an ESP in front of a Sun-like star. 
Shown are results for two planet radii, Rp = 1 Rj^p (black) and Rp = 2 Rj-^p (blue). The solid 
lines correspond to central transits {i = 90°), the dashed lines to an inclination of i = 87°, 
and the dashed lines to an inclination of i = 84°. The planet orbital radius is 0.04 AU. 
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Fig. 4. — Occultation polarization at 4600A for the transit of an ESP in front of a Sun-like 
star. Shown are results for two planet radii, Rp = 1 Rj^p (black) and Rp = 2 Rj^p (blue). 
The solid lines correspond to central transits {i = 90°), the dashed lines to an inclination of 
i = 87°, and the dashed lines to an inclination of i = 84°. The planet orbital radius is 0.04 
AU. The horizontal line indicates the best polarization sensitivity for current broad-band 
polarimeters. 
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Fig. 5. — Same as Fig. 4 for the resonant line Ca I 4227A. 



Table 1. Limb Darkening Coefficients 



Spectral Type 


ai 


02 


as 


04 




Teff (K) 


G2 (sun @ 4227A) 


0.q450 


-0.256 


1.286 


-0.584 


1 


5780 


G2 (sun @ 4600A) 


0.368 


0.265 


0.588 


-0.356 


1 


5780 


K5 (l^-band) 


0.419 


0.704 


-0.220 


-0.035 


0.72 


3500 


M2 (l^-band) 


0.396 


0.907 


-0.402 


0.086 


0.5 


3000 


T (l^-band)'^ 


0.396 


0.907 


-0.402 


0.086 


0.2 


3000 



^To our knowledge there is no published data about limb darkening of T 
dwarfs. For this reason, we used the limb darkening coefficients of the M2 
dwarf, as an approximation. 
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Fig. 6. — Limb polarization for different values of k. Black: k = 0; blue: k = 1; green: 
k = 5; and red: k = 50. The dotted line represents the curve for the solar limb polarization 
at 4600A and the dashed line represents the limb polarization for a M giant (Harrington 
1969). 



- 20 - 



0.025 




t (h) 



Fig. 7. — Occultation polarization efficiency for the central transit of a 1 fZjup planet accross 
a M2 dwarf. Black: A; = 0; blue: k = 1\ green: k = 5; and red: k = 50. For Pi ^ 1%, the 
OP should be measurable if e > 0.001. 




Fig. 8. — Maximum polarization efficiency vs. planet radius for a M2 dwarf. Black: k = 0; 
blue: k = 1; green: k = 5; and red: k = 50. The dotted lines represent the analytical 
approximation of Eq. 11, with the appropriate values of w, k and Oj. 
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10. — Same as Fig. 8 for a K5 dwarf. 
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Fig. 11. — Same as Fig. 7 for a T dwarf. 
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12. — Same as Fig. 8 for a T dwarf. 
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Table 2. Observed Limb Polarization Coefficients for the Sun 



Wavelength (A) 


Pi 


k 


Reference 


4227 


0.165 


50 


Bianda ct al. (1999) 


4600 


8.7 10~3 


70 


Faurobcrt ct al. (2001) 


5160 


3.3 10-3 


59 


Faurobert & Arnaud (2003) 
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Table 3. Maximum OP for Representative Models of Late Type Stars, using Pi = 10% 



Spectral Type -Rp = 1 -REarth -Rp = 1 ^Earth -Rp = 2 R^^p -Rp = 2 Ry 
k = 5Q k = Q k = 50 k = 



K5 7.5 X 10-^ 1.1 X 10-5 1.0 X 10"^ 2.6 x 10"=^ 

M2 1.7 X 10-^ 2.4 X 10-5 1.7 x 10"^ 4.6 x 10"=^ 

T 9.0 X 10-6 1.6 X 10-^ 3.9 x 10-^ 9.2 x 10"=^ 



